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Periodic solution parameters, in chaotic dynamical systems, form periodic windows with charac- 
teristic distribution in two-parameter spaces. Recently, general properties of this organization have 
been reported, but a theoretical explanation for that remains unknown. Here, for the first time we 
associate the distribution of these periodic windows with scaling laws based in fundamental dynamic 
properties. For the Rossler system, we present a new scenery of periodic windows composed by mul- 
tiple spirals, continuously connected, converging to different points along of a homoclinic bifurcation 
set. We show that the bi-dimensional distribution of these periodic windows unexpectedly follows 
scales given by the one-parameter homoclinic theory. Our result is a strong evidence that, close to 
homoclinic bifurcations, periodic windows are aligned in the two-parameter space. 
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I. INTRODUCTION 



For several smooth nonlinear maps and differential 
equations, stable periodic orbits and their dependence 
on the system control parameters are well known. The 
existence of these orbits can be properly visualized on a 
bi-dimensional parameter space, where generally we find 
periodic windows, i.e., continuous sets of parameters, em- 
bedded into chaotic regions, for which periodical orbits 
exist [U [5] . There is an intricate periodic window, quite 
general in dynamical systems, whose local nature was 
explained in [Tl [3] but the global features, despite the 
great number of studies [4lll6j. remain not well under- 
stood. For local and global features we mean typically 
qualities of an isolated and a multiple periodic windows, 
respectively. Recently, it has been found that these peri- 
odic windows, baptized shrimp [2] , are continuously con- 
nected along spirals emerging from in a homoclinic bi- 
furcation point [T71 HH] (set of parameters for which a 
bi-asymptotic curve, the homoclinic orbit, converges to a 
saddle- focus equilibrium point). These spiral structures 
were verified experimentally in [121 HO] and are also ob- 
served in d mills • 

However, in these researches the distribution of 
shrimps in the parameter space have not yet been clearly 
associated to any fundamental dynamical property. To 
accomplish this, we investigate the relation between the 
shrimps and the homoclinic curves in the parameter 
space [231 121] [30] . For the Rossler system, we present a 
new and remarkable two-parameter space scenery where 
from each shrimp emerge infinity spirals with focus in dis- 
crete points along of a homoclinic bifurcation curve (con- 
tinuous parameter sets for which homoclinic orbits exist). 
Each spiral is composed by a shrimp family, i.e., infinite 
shrimps continuously connected in a spiral sequence. We 
show that, even the shrimps are a codimension-two phe- 
nomena (two parameters are necessary to obtain it), they 
are accumulating at the spiral focus following scaling laws 



predicted by the one-parameter space homoclinic theory 
[25] . We also show that the reported period adding cas- 
cades observed in shrimps accumulations [HI [12] is a con- 
sequence of the spiral periodic windows approach to the 
homoclinic bifurcation point. 

Recently was published a work [20] about scaling laws 
in a electronic homoclinic system where the authors as- 
sociate scales of tangent bifurcation with the shrimp dis- 
tribution. Here we discuss this issue in details and show 
that, in the two-parameter space, the scales measured 
in shrimps correspond to the distance between crosses of 
superstable periodic curves. Furthermore, we call the at- 
tention to the necessity of a rigorous prove of these scales 
in shrimp distributions which was done in [26j . 

In section JTTjwe present the spiral scenery of periodic 
windows, in |HI| we show the scaling laws concerning its 



distribution and periodicity, and in section IV we present 
the conclusions. 



II. SHRIMP DISTRIBUTIONS 



We consider the homoclinic Rossler system given by 



X = —y — z 
y = x + ay 
z = bx — cz 



(1) 



where we fix & = 0.3 and analyze the parameter space c' x 
a', where a — a' sm{6) + c' cos{9) + gq and c = a' cos{9) — 
c'sm{e) + Co with ao = 0.3301, co = 4.9305, and 9 = 
88.8°. For the region investigated (Fig. [T]) the origin 
phase space (Pq) is a saddle- focus and the eigenvalues 
of Eqs. (IT]) Jacobian matrix evaluated at Pq are \i 2 — 
pzt iuj and A3 = A, where A, p, and w are M*. 

Periodic and chaotic asymptotic solutions of Eqs. ([T]) 
are determined numerically by evaluating the largest 
nonzero Lyapunov exponent I. In Fig. [T] from black 
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FIG. 1; (Color online) (a) Parameter space c' x a . The phase space has an chaotic attractor for parameters in yellow and 
periodic for parameters in gray. The sets A2 to are shrimps connected in a spiral way converging asymptotically to the 
onset of H\ homoclinic curve (dashed line). B, C, D, and E are shrimp families emerging by flip bifurcations in A2 and A3. 
The blue square regions are magnified in (b), (c), and (d). The black dots are the spiral focus and the blue dots indicate the 
shrimp positions. 
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FIG. 2: (Color online) Schematic scenery. Spiral in black corresponds to family A. In red are spirals emerging from shrimps 
A2 and A3. The dashed line is the Hi homoclinic bifurcation. The dots along Hi are the spiral focus points while the points 
in the spirals represent the crosses of superstables of A2, A3, Bo, Co, Do, and Eq. 



to white {I < 0) the periodic orbits increase their stabil- they are superstable {I achieve the most negative value), 
ity: in black these orbits bifurcate (/ = 0) and in white Otherwise, from black to yellow {I > 0) the behavior is 
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asymptotically chaotic. The structure labeled A3 is a 
particular region of the parameter space where a typical 
periodic window, the shrimp, can be visualized. The big 
shadow region, from which extends four narrow antennae, 
is its central body and corresponds to the fundamental pe- 
riodic window where is the fundamental periodic orbit of 
the shrimp |15) . The central body is bordered by tangent 
bifurcations, where we find chaotic regions, and by flip 
bifurcations, where we find sequences of doubled periodic 
regions with similar shape to the central body. The white 
lines form the shrimp skeleton (parameters set where the 
periodic regions are around) and represent superstable 
periodic orbits with behavior strongly attractive. Note 
that, in the central body, there is a cross between two 
remarkable superstable curves (blue dots in Fig. [T]). We 
consider this point as a shrimp position and localize it 
identifying these intersections. In this sense, shrimps are 
codimensio-two structures. 

To obtain the parameter sets of the homoclinic bifur- 
cation Hi associated with the saddle- focus point Pq, we 
determine numerically the parameters for which the sta- 
ble and unstable manifolds of Pq merge constituting a 
homoclinic orbit (To obtain homoclinic orbits in piece- 
wise systems, see Ref. f27]). The labels A, B, C, D, 
and E, in Fig. [TJ indicate different shrimp families asso- 
ciated with Hi. Note that, in each family, shrimps are 
continuously connected along a spiral sequence around 
a homoclinic bifurcation point [Typically, two antennae 
of each shrimp central body is connected with the pre- 
vious and the next shrimp and the other two antennae 
connect the shrimp to the divergence region (a' > 0.2)]. 
The index i, in the label families, indicates the element 
of these sequence. The first shrimp of the sequence is 
indexed hy i — I, the second by i = 2, and so on. The 
index i = indicates the structure that connect different 
families. Thus, Bq connect family A with family B. Fur- 
thermore, excluding the central body of A2, for example, 
Bq has the same features that characterize family B. 

We observe that each shrimp generates many others 
spiral sequences by flip bifurcations. The central bodies 
of A2 and ^3, in Fig. [ija), bifurcate in Bq, Cq, Dq, and 
Eq, which connect the family A with the families B, C, 
D, and E. The family B converges to the region close the 
point (a' = 0.0228, c' = 0.0550) in an anti-clockwise di- 
rection and the family C converges to the region close the 
point (a' ~ 0.0333, c' = 0.0800) in a clockwise direction. 
Both focus seems to be along the Hi homoclinic bifur- 
cation curve [see magnifications in Figs, [ijb) and (c)]. 
The families A, D and E converge to the Hi onset point 
(0,0) [see Figs. fT^a) and (d)]. The schematic scenery 
is shown in Fig. 2j As far we could check, the same 
is observed in any element of any family suggesting the 
existence of a fractal structure of spiral self-replications 
converging asymptotically to the homoclinic bifurcation 
curve. Moreover, we identify shrimps concentrated in two 
sets, between the two shrimp antennae that converge to 
the divergence region (as B and C families), and between 
two shrimps, i and i -|- 2, of the same family (as D and 



E families). 

With respect the trajectory behavior, following contin- 
uously a shrimp spiral sequence, we observe that the orbit 
time-period grows smoothly tending to infinity close the 
homoclinic orbit. So that, the fundamental periodic orbit 
add one cycle from i-shrimp to the (i + 2)-shrimp form- 
ing a period adding cascade accumulating into the Hi 
curve, as shown in Figs. [3] (a)-(d). The period adding 
can be identified considering a flat Poincare section in 
X — with y < close to the plan xy (represented by 
the dashed lines). The orbits cross the Poincare section 
two times in (a), three times in (b), and four times in 
(c). The homoclinic orbit Hi has infinite cycles in (d). 
Note that here we considered Aq as a shrimp as discussed 
before. 



III. SCALING LAWS FOR SHRIMP 
DISTRIBUTIONS 

Next we analyse our numerical results presented in 
the last section from the homoclinic theory described 
by Shilnikov theorem. We focus in the characterization 
of the two-parameter shrimp structures from these one- 
parameter theory. 

Shilnikov theorem can be applied to systems which 
saddle-focus equilibrium point have homoclinic orbits (F) 
solutions and, in its normal form, can be represented by 



x' = px' - ujy' + 0{x', y', z') 
y' = ujx + py + P(x ,y ,z ) 
z' = ~\z' + Q{x',y',z'), 



(2) 



where O, P, and Q are analytic functions with 0(a?') = 
P{x') = Q{x') = 0, for the equilibrium saddle-focus point 
x' — 0, and A, p, and uj are M^. 

If the saddle-focus form F with saddle index f < 1/2 
{ly = p/A), the Shilnikov Theorem shows that, close to 
F, in a one-parameter space, there are infinite countable 
sets of periodic and homoclinic bifurcations accumulating 
into F, namely: (I) Stable periodic solutions emerge dis- 
tributed as Sp = limi_^oo(Mi-(-i/Mi) = exp — (ttp/w) where 
Pi and pi+i are two consecutive tangent bifurcation pa- 
rameters [H [25l [28| • The period difference between two 
consecutive periodic orbits is AT = \mii_yao{Ti^i — Ti) — 
■n /lo, where Ti is the period of the «-periodic orbit; (II) In- 
finite classes of homoclinic orbit, characterized by similar 
orbits, called secondary homoclinic orbits [double-pulse 
(r2), triple-pulse (F3), and n-pulse (r„)] are distributed 
following 5'r„ = limi_j.oo(Mi+i/MD = exp — (ttA/cli), where 
p'^ and p'ij^i are two consecutive homoclinic bifurcation 
parameters of F„ |29| . The limit i = 00 is related with 
the primary single-pulse (Fi). 

To associate the shrimps with this theory, we verify 
these scaling laws in our numerical simulation. For the 
family A, the homoclinic orbit parameter is (a^,c^) = 
(0, 0) for which the eigenvalues calculated at the fixed 
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FIG. 3: Period adding of A shrimp family. The periodic trajectories cross the Poincare section (dashed Hnes) two, three, four, 
and infinity times from (a) to (d), respectively. 



point Po are Ai,2 = 0.1354 ± iO.9867 and A3 = -4.8713. 
Thus, the scaling laws are Sp = 0.6498, AT = 3.1840, 
and S'r ~ 2.10"^. To measure numerically the pe- 
riodic orbit scaling Sp, we define Si — di+i/di, with 



+ (5^+1 Where (a^,c^) 



intersection point between two superstable curves as de- 
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FIG. 4: (Color online) (a) and (b) show, respectively, measure 
of distribution and difference time-period scalings of shrimps. 
The dashed lines are the theoretical values, Sp = 0.6498 and 
AT = 3.1840, in the limit i -s> 00, of the family A. The 
symbols square, triangle, and circle are respect families A, B, 
and C. 



fined before the i-shrimp position. The measured Si for 
the shrimp family A converges quickly to the theoreti- 
cal scaling parameter Sp as shown in Fig. |4][a) (curve 
with square symbol). Figure |4]^b) shows that the mea- 
sured time-period differences AT^ = Ti+i — Ti also con- 
verge quickly to the theoretical scaling AT. Thus, in 
the limit i — !■ 00, the period difference between Ai and 
Ai^2 is 27r/w = 6.3681. This interval time corresponds 
to one periodic orbit revolution close to the unstable 
manifold (plan xy) and explains the period adding phe- 
nomenon in homoclinic systems. It suggests that sim- 
ilar mechanism can explain the period adding in other 
systems. For the family B, we consider {a'^,c'^) = 
(0.0228,0.0550). The eigenvalues calculated at Po are 
Ai,2 = 0.1524 ± i0.9831 and A3 = -4.7799. Thus, the 
scafing laws are Sp = 0.6144, AT = 3.1956, and Sr w 
2.10~^, similarly to the family A. And, for the family 
C, ia'^,c'^) = (0.0333,0.0800), A1.2 = 0.1522 ± i0.9831 
and A3 = -4.7910 with Sp = 0.6148, AT = 3.1955, and 
Sr « 2.10"'^. The scalings S", and AT, of B and C fami- 
lies are in agreement to the theoretical estimated Sp and 
AT as shown in Fig. |4f^a) and (b) (triangle and circle 
symbols for B and C families, respectively). Although 
the homoclinic scaling law parameters is very low for be- 
ing verified numerically, we have observed the existence 
of many different bifurcations i7„ close to the primary 
Hi homoclinic curve parameter. The scalings measured 
of families D and E have values very similar to the family 
A and it was omitted in Fig. |4] 
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We emphasize that it is not expected that the scahng 
laws here considered match with the shrimps scahngs, 
since shrimps are a codimension-two phenomena. Our re- 
sults suggest strongly that shrimps are distributed along 
lines in regions close to the homoclinic orbit bifurcation 
[See Figs.[l](b)-(d)]. 

IV. CONCLUSIONS 

We used the knowledge of homoclinic orbits distribu- 
tion in the bi-dimensional parameter space to explain the 
distribution of periodic windows in this space. We found 
infinite periodic structures with spiral shape distributed 
along the homoclinic bifurcation curve Hi. Each spi- 
ral has two extremities: one is a point in Hi, which is 
its focus, and the other is a shrimp. In the other side, 
from each shrimp derive infinite spirals with different fo- 
cus along Hi. Thus, the spiral are composed by infinite 
shrimps, which compose infinite spirals which have infi- 
nite shrimps which compose infinite spirals and son on. 
This characterize the fractality of this scenery of spirals. 

It is worth to mention that, excluding family A, each 
spiral is composed by a family of shrimps that, close the 
Hi curve, the shape of its orbit in the phase space ap- 



proaches the shape of an secondary homoclinic orbit r„ 
(also called subsidiary homoclinic orbit). It suggest that 
the real focus of each spiral is a point close Hi where is 
formed a r„ orbit. As discussed before, the classes of ho- 
moclinic orbit are organized close Hi following the scale 
5r„ • Thus we argue that the spiral distribution in the 
parameter space should follows the scale Sr„ . 

We have also identified properties of the homoclinic 
theory in the shrimps organization. We notice that this 
result is not expected since the homoclinic scalings are 
valid just in a one-parameter space and shrimps are two- 
parameter structures. It indicates that shrimps, close 
to the homoclinic bifurcation, are organized along a line 
that intersects Hi . The analytical prove of this is merit 
of investigation [26 . 
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